The cosmological gravitational wave background from primordial density 

perturbations 
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We discuss the gravitational wave background generated by primordial density perturbations 
evolving during the radiation era. At second-order in a perturbative expansion, density fluctua- 
tions produce gravitational waves. We calculate the power spectra of gravitational waves from this 
mechanism, and show that, in principle, future gravitational wave detectors could be used to con- 
strain the primordial power spectrum on scales vastly different from those currently being probed 
by large-scale structure. As examples we compute the gravitational wave background generated by 
both a power-law spectrum on all scales, and a delta-function power spectrum on a single scale. 



I. INTRODUCTION 

Gravitational waves are an inevitable, yet still elusive, 
consequence of Einstein's theory of general relativity that 
will be tested and, we hope, revealed by upcoming experi- 
ments. Linear perturbations about a cosmological metric 
include transverse and trace-free modes (tensor modes) 
that propagate independently of conventional matter per- 
turbations at first order - i.e., gravitational waves. Dur- 
ing an inflationary expansion in the very early universe 
large scale (super-Hubble scale) tensor modes will be gen- 
erated from initial quantum fluctuations on small scales 
EL S Hi • But the amplitude depends on the energy scale 
of inflation and may be unobservably small if inflation 
occurs much below the GUT scale 

On the other hand primordial density perturbations, 
and their associated scalar metric perturbations, do exist 
and these will inevitably generate a cosmological back- 
ground of gravitational waves at second order through 
mode coupling 0, @, 0, & 1 E3- 

Given the detailed 
information we now have about the primordial density 
perturbations on a range of cosmological scales, it is now 
timely to consider the amplitude and distribution of ten- 
sor (and vector) modes that will be generated at second 
order. 

The observed Gaussian distribution of the primordial 
density perturbations will generate second-order modes, 
which will have a x 2_ distribution, unlike any first-order 
gravitational waves from inflation. Given the observed 
primordial power spectrum of scalar perturbations on 
large scales, of order 10~ 9 11], one would expect the 
power spectrum of second-order metric perturbations to 
be of order 10 -18 in the radiation dominated era. But 
this naive expectation needs to be tested against a full 
second-order calculation. To our knowledge this has not 
previously been done. Recently the effect on the cosmic 
microwave background (CMB) of second-order gravita- 
tional waves generated on very large scales was inves- 
tigated by Mollerach et al 12] but these modes enter 
the Hubble scale after matter-radiation equality. We will 



consider much smaller scale modes that enter the Hub- 
ble scale during the primordial radiation-dominated era 
and may be relevant for direct detection by gravitational 
wave experiments. 

Moreover, the generation of gravitational waves from 
primordial density perturbations on smaller scales, not 
directly probed by astronomical observations, could be 
used in the same way that primordial black hole forma- 
tion has previo usly been used to constrain overdensitics 
on these scal es fl3f . Recently Easther and Lim [l4| (see 
also Ref.pjJ lid. Il7|) have suggested that large density 
inhomogeneities on sub-Hubble scales during preheating 
at the end of inflation could generate a gravitational wave 
background that might be detectable by future gravita- 
tional wave experiments such as Advanced LIGO (l8j . 

In this paper we present the second-order evolution 
equation for gravitational waves generated from terms 
quadratic in the first-order matter and metric perturba- 
tions. In section |TT] we present the field equations for 
the background Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metric and first-order perturbations, giving the 
standard solutions in terms of the Bardeen metric poten- 
tial during the primordial radiation-dominated era. In 
section IIIII we present the second-order evolution equa- 
tion for gravitational waves driven by a source which is 
quadratic in the first-order scalar perturbations. We use 
the first-order constraint equations to eliminate the mat- 
ter perturbations and write this evolution equation solely 
in terms of the Bardeen potentials and their derivatives. 
Finally we present solutions to the second-order gravita- 
tional wave equation using a Green function method. We 
calculate the power spectrum of the gravitational waves 
generated first by a delta-function spectrum of density 
perturbations at a particular wavelength, and then by 
a power-law spectrum for the primordial density pertur- 
bations. We conclude by discussing the expected ampli- 
tudes of the gravitational waves produced and compare 
with sensitivities of current and future detectors in sec- 
tion nvi 
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II. DENSITY FLUCTUATIONS IN THE 
RADIATION ERA 

We will consider perturbations up to second order 
about a Robertson- Walker background. We decompose 
the metric as 1331 



laf} = 9af} + Sg a p + S 2 g a/ 3. 



(1) 



As we are considering gravitational waves sourced by 
first-order scalar perturbations, 5g a p has purely scalar 
degrees of freedom, while 5 2 g a p has, in general, scalar, 
vector and tensor modes induced by 8g a p. However, we 
are only investigating second-order tensor modes so will 
project out any scalar or vector modes at second order; 
we may therefore consider 5 2 g a p as having purely ten- 
sor degrees of freedom. Choosing a longitudinal gauge at 
first order, we write the metric as 



5oo 



-a 2 [1 + 2$] , 



(l-2*) 7 y 



9oz 

I — hi 4 
2 3 



0, 



(2) 



Here, $ is the first-order Bardeen potential, and hij is 
the second-order tensor mode. Note that 7 J /lyifc = 
and j l ° hij = 0. 

We assume a spatially flat geometry and a pure ra- 
diation background. The scale factor, Hubble rate and 

energy density evolve as a — do H = a ^ = V 1 

and, p oc fj^ 4 , in terms of conformal time rj. 

At first-order, assuming no anisotropic pressure, the 
Bardeen potential, for a comoving wavenumber k 7 satis- 
fies 



+ m (1 + c 2 ) $' + c 2 s k 2 <S> = 0, 



(3) 



where the speed of sound is c 2 — 1/3, and a prime denotes 
a derivative with respect to conformal time. The general 
solution to this is 



$(fc,77) = 



[ki]) 3 

B(k) 

{krj) 3 



krj 

71 

krj 

7i ! 



cos 



krj 

71 

krj 

71 



krj 

71 

krj 

71 



(4) 



This will act as a source for the GWs at second order. 
At early times, krj — > we see that 



$(fc,7?) 



m B(k) 

9V3 (kr/) 3 



(5) 



the second term is the decaying mode which we shall 
neglect hereafter. 

Assuming that the fluctuations are Gaussian, we may 
write A(k) — A(k)E(k) where the E are Gaussian ran- 
dom variables of unit variance which have the property 

(E*(kx)E(k 2 )) = 5 3 (k 1 -k 2 ). (6) 
The power spectrum for the scalar perturbation can then 
be defined as 

($*(fci)$(fca)) = ^-5(fc! - fea) (7) 
implying that at early times the power spectrum becomes 

A- 3 



486tt 2 



(8) 



The Bardeen potential can be related to the comoving 
curvature perturbation at early times, giving 



A{kf 



21671-' 

k 3 



(9) 



where A 2 ^ is primordial power spectrum for the curva- 
ture perturbation 1Z. Current observations show A\ ~ 
2 x 10~ 9 at a scale kcuB = 0.002Mpc _1 , and is almost 
independent of wavenumber on these scales [ll[ . 



III. THE INDUCED GRAVITATIONAL WAVES 

We now consider the evolution equations for the sec- 
ond order tensor perturbations, hij, sourced by the scalar 
density perturbations discussed above. 

We will write the Fourier transform of ha as 



hij(x,r)) 



(2tt) 3 / 2 
+ h(k,rj)qij(k)] , 



d?ke lk x [h(k,rj) qij (k) 



(10) 



where the two polarization tensors and q~ij are ex- 
pressed in terms of the orthonormal basis vectors e and 
e orthogonal to k, 



qij(k) 
qijik) 



1 

71 
1 

71 



\e l {k)ej{k) - e l {k)ej{k)} , 
[ ei (fe)e,-(fe)+ei(fe)ej(fe)]. (11) 



Thus to extract the transverse, trace-free part of any tensor we project with the operator Tij lm , defined through its 
action on a two-index tensor 



7"' lm c 



d 3 k' 



(2tt) 3 / 2 



Am 



d 3 x < ^Ke^-*')^^') 



(27T) 3 / 2 



d 3 k> 



(2tt) 3 / 2 



a , q lm (k') 



d 3 x' 



(2tt) 3 / 2 



ih • (a 



] S lm {x'). 
(12) 
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The evolution equation is calculated by expanding the Einstein field equation up to second order including quadratic 
terms in the first-order scalar modes. The only contributions from pure second order matter fluctuations would come 
from a transverse and traceless contribution to the anisotropic stress which we ignore in this analysis. Thus, calculating 
the transverse, trace-free spatial part of the field equations yields 0, [2(| 



4 + 2Hh' lJ - V 2 hj = -4%/ m S lr , 



The source term is given by [l9l |20T | 



2 2 



n 2 ^ b + 2n^ b + 



(13) 



(14) 



A pipe denotes the spatial covariant derivative, and k 2 — 8ttG. This expression is consistent with the second-order 
Ricci tensor for scalar perturbations about an FRW background calculated in Ref. [2l| . 

In Fourier space we then find the amplitude of the tensor mode, for either polarisation, obeys the evolution equation 



h"(k, n) + -h'(k, 77) + k 2 h(k, rj) = S(k, rj), 



where the source term is given by 



S(k,fj) 



d i kk i k j 12$(fc - k,fj)$(k,f)) 



~2 

rf 



f)<S>{k - k, rj) + (k - fc, 77) 



*'(M) 



(2tt) 3 / 2 

The particular solution for the gravitational waves is then given by an integral over the Greens function 



h(k,rj) = 



1 



a{rj) 



Gk{r],fi)a{fj)S{k,fi)dfi, 



'la 



where 



(15) 



(16) 



(17) 



(18) 



Gk{"n,fj) = —57- sm(kri) cos(kfi) — cos(krj) sin(fcry) 

The power spectrum of the induced GW is defined in the usual manner, 

2tt 2 

(h(ki,r))h(k 2 ,ri)) = ^r^( fe i - fc 2) V h {k, rj). (19) 
Substituting for h(k,Tj) and using Wick's theorem, and using spherical coordinates in Fourier space we find that 



V h {k,rj) 



1 



8n 4 a(r)) 2 
The variable u is given by 



I/O 



df] 2 I dfjx a (771)0(172) 



dk 



du /c 2 G /c (t7, 771)6^(77, 77 2 )F(fc,fc,w; 771, 7)2). (20) 



1 + (k/k) 2 - 2(k/k) cose, 



(21) 



where 9 is the angle between the modes k and fc, and the limits of integration u± correspond to the angles 8 — 0,tt 
respectively. The integrand F in Eq. (12"U1) is found to be 



F{k,k,u;rix,fi2) = uk (2kk) 2 - {{ukf - k 2 - fc 2 ) 2 



|3$(77fc, 77i)$(fc, 771) + [277i$(7tfc, 771) + fj 2 $'(uk, 7)1)] $'(fc, 77i)| 
|3$(itfc, 7) 2 )$(fc, 772) + 77 2 $'(wfc, 77 2 )$'(fc, 772) + 772 <&{uk, 7)2) + $(fc, 77 2 )$'(7ifc, 772) V 



22) 



where we distinguish between the Gaussian random variable $(fc, 77) in Eq. (| 16[) and its amplitude <&(&, 77) in Eq. ([22 
which we take to be isotropic. We introduce the following dimensionless variables 



k' 



x = kr\ 



(23) 
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and substitute for the first-order solution for <f> from Eq. (|4|), as well as the Greens function, into our GW power 
spectrum. After some simplification and substitution from the above formulas we have 



2(216) 2 
ttV 



dv 



\v+l\ 



du 



1 



(uvy 



Uv 



2 (u 2 - v 2 - l) 2 ] 2 V^(uk)V^(vk) 



sin(a;) / dx\ I\{x\) — cos(x) I dx\T?.{x\) sin(a;) / 

J Xn J Xn J L J X 



sin(a;) / dx^ Ii(x2) — cos(x) I dxil^x^) 



(24) 



where we have defined the four functions 



M 3 



x ' 



(25) 



m— 1 n— 1 

The coefficients a n ,<p n and Mij lm in this expression are dependent on u and v but not x, and may be found in the 
appendix. 

Evaluating the integrals in the rhs of Eq. (|24p for various input scalar power spectra will then tell us the power in 
each GW mode. This is not particularly simple, so we start by analytically expanding the integrals over the functions 
Xj by parts up to Si and Ci functions [22l |: 



Xj(u, v, x, xq 



rx 5 8 

)= / dxljix) = E E M » 
J Xn 1 



m— 1 n— 1 
(m-1) 

(to - 1)! 



(m-fc-2)! fc . / . , (fc + 2) 
2 y • i — ; ttt^-o;!; sin [ a n x 1 



fe=i 



(to - 1)! 



7T I x( 1 + fc -™) 



da; — sin a n x + 0„ 

x 



(m + 1) 



(26) 



The remaining two integrals over Fourier space can now be done numerically once power spectra for the scalar modes 
are chosen. We shall only consider modes which start their evolution well outside the Hubble radius, and hereafter 
set xq = 0. 



A. Gravitational wave generation by a single scalar 

mode 

Second-order gravitational waves potentially provide 
a method by which we could detect a particular scalar 
mode with excessive power, compared to the roughly 
scale-invariant average we observe on large scales today. 
As a precursor to a power-law power spectrum for the 
scalar modes, we can investigate how the system reacts 
to power being put in at one particular scale. That is, we 
put power in at a single wavelength (i.e., a single comov- 
ing scale) which is described by an isotropic distribution 
of wavevectors (i.e., at all possible angles). This would 
be useful when considering preheating for example which 
can result in features (large power over narrow range of 
scales) in the power spectrum of the input scalar pertur- 
bations 0Q3]. 

We will therefore consider the case of the delta- 
function power spectrum, as an idealised limit of a spike 
in the power spectrum, or of power being introduced on 
narrow range of scales above the roughly scale-invariant 
primordial spectrum observed. The specific form we 
choose is: 

P 9 (k) = ^A 2 A 2 n (k CMB ) 5(k - fc m ), (27) 
where A is the amplitude at a single wavenumber, ki n , 



relative to the observed amplitude of the primordial 
power spectrum, A^kcMB), at wavenumber kcMB 3> 
k 

The power spectrum of the gravitational waves pro- 
duced by this scalar mode is then, from Eq. I] 24 p. 

2(2\6) 2 r°° /•l t '+ 1 l 

P h (k,n) = / dv' / du> 

* V JO J\v-1\ 

V<s,{u'k) V$(v'k) T 5 {u',v',x), (28) 

where 

F S {u,v,x) = -±-,[Av 2 -(u 2 -v 2 -I) 2 ] 2 
(uv) 8 L J 

x [sin(x)Ai — cos(x)A 2 ] [sin(a;)A"3 — cos(a;)X4] .(29) 

We may now evaluate the v! and v' integrals for a delta- 
function power spectrum, giving 

Vh(k,v) = ^tt(^ R ) 4 ^(v,4 (30) 
provided that the following condition holds: 

»=»=¥4 < 3i > 

The inequality arises from the fact that gravitational 
waves cannot be excited with more than twice the mo- 
mentum of the density perturbations, fej„. 
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20 40 60 80 100 



FIG. 1: The function J-s(x) for different values of ki n /k = 
u = v, with xq = 0. 



Although .7"a(t7, u, x) is a very complicated function to 
write down, its properties are relatively straightforward 
to understand, as we can see from Figs.Q]-[D 

In Fig. Q] we see the basic behaviour for small x. When 
xq = 0, we may interpret x as time for a fixed wavelength 
or wavenumber at fixed time. The function therefore is 
zero at early times and large scales, has interesting oscil- 
latory behaviour on scales of order the horizon size, when 
krj = 1, and oscillates at constant amplitude at small 
scales and/or late times - i.e., when modes are well inside 
the horizon. The asymptotic behaviour is clearly seen in 
Fig. [2j where we are looking at modes with fc = fcj n , 
generated when there is a 60 degree angle (9 = 60°) be- 
tween the input and output mode (where 9 is given by 
Eg . |2"T1 and Eq. l23[) . We have power-law growth for small 
x, and constant amplitude for x — > oo, with interesting 
oscillations on scales of the order of the Hubble scale. In 
Fig. [3] we show the same but now with u = v = vo/2. 
We see that the amplitude continues to grow on small 
scales or late times, corresponding to resonant amplifica- 
tion of modes generated at an angle cos# = l/y/3 to the 
incoming mode. 

In Fig. Q] we show the envelope of J-g at fixed large 
x (where we have chosen x = 10 6 ) which shows the 




-1 1 2 3 4 5 6 



log l0 x 

FIG. 2: The function Ts[x) from early times (super-Hubble 
scales) to late times (small scales) for k = kin (i.e. u = v = 1), 
xo = 0. For large x the function oscillates with a constant 
amplitude. 




-6 



2 4 6 8 10 12 14 

log 10 x 

FIG. 3: The function Ts(x) from early times (super-Hubble 
scales) to late times (small scales), for the resonant case of 
ki n /k = u = v = VS/2. For large x the power in the gener- 
ated GW continues to grow logarithmically. 



shape of the power spectrum one could observe at late 
times. There is a power-law tail on large scales, k <C ki n 
which is proportional to (fc/fcin) 3 as one would expect 
for perturbations generated from much smaller scale 
modes and hence uncorrelated on larger scales. On 
small scales there is a sharp high-frequency cutoff at 
k = 2ki n . We also see the resonance as a sharp spike in 
the power spectrum, of approximate width 1/x and am- 
plitude 0.05(log 10 x)(log 10 0.02a;). Finally we note there 
is zero power at u — v = y/3/2, corresponding to an angle 
cos 9 = l/v6 between the input and generated modes. 
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log w k/k in 

FIG. 4: The power at late times, x — 10 20 , for all scales, 
showing a resonance at u = V3/2 and a power-law tail (due 
to modes on larger scales being uncorrelated) . Detail of the 
resonance is shown in the inset revealing a wobbly structure. 

B. Power law scalar modes 

While the preceding subsection gave insight into many 
aspects of the generation mechanism, we would also like 
to know the GW generated from a nearly scale-invariant 
spectrum of density fluctuations. To investigate this, we 
assume that the input power spectrum is: 

*<*->- 5** G^r'- < 32 > 

where the index n s tells us the tilt of the spectrum rel- 
ative to scale-invariance, n s = 1, and kcMB is a pivot 
scale for the power spectrum [TT| . The power spectrum 
of the generated gravitational waves is then given, from 
Eq.OUby 

r h (k,v) = 2(21 f 2 A ^ (^-Y^ r n .{x). (33) 

K V \KGMB J 

The function J-„ a is defined by 

JF„ s (x) = / dv du (to)"* -1 !Fs(u,v,x). (34) 

Jo J\v-1\ 

We show a plot of this function for a scale-invariant spec- 
trum, n s — 1 in Fig. [5l 

We can think of T ns (x) as a function of wavenumber 
k = x/rj at a specific time rj, in which case the amplitude 
°f 3~n s peaks on scales just inside the Hubble radius, and 
becomes scale-invariant on smaller scales. The series of 
oscillations for x > 1 is analogous to the acoustic peaks 
in the CMB spectrum, although one might have expected 
this to be less pronounced at second order where there 




1 10 100 

X 



FIG. 5: Tn s (x) for a scale invariant input power spectrum. 
For this case the power spectrum of gravitational waves is also 
scale invariant owing to the scale invariant conversion factor 
between T na and Vk given by Eq. (p3)l . 

is an integration over modes. It suggests that a narrow 
range of lst-order modes with well-defined phase make 
the dominant contribution to the second-order tensors 
on a given scale. 

On the other hand, Fig.[5]also represents the evolution 
of a single mode over time. Thus, power is continually 
added to the GW until it is well inside the Hubble radius, 
before oscillating at almost constant amplitude at late 
times, i.e., a freely propagating GW. 

How does the tilt of the scalar modes affect the gen- 
erated GW? As far as T ns is concerned there is only a 
small amplification change, shown in Fig. [6] for large x. 
There is also an increase in the power at small scales due 
to the factor {k/kcMB) 2 ^"^ 1 ^ in Eq. (33]if we consider a 
blue spectrum (n s > 1). 




0.8 0.9 1 1.1 1.2 



n s 

FIG. 6: J- ns for x — > oo, xo = 0, vs. the tilt n s of the scalar 
power spectrum. We have normalised the vertical scale to 
JF„„ = i = 8.3 x 1(T 3 . 

Finally, we note that for non-zero xo there will be a 
downward break for large k in the scale invariant F ng , 
which corresponds to kr]Q ~ 1 - i.e., modes which are 
inside the horizon at the start of the interaction period. 
However, if we start our integration at the GUT scale this 
will only be relevant for very high frequencies (^ 10 8 Hz), 
so we shall not pursue this further here. 
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IV. DISCUSSION 

We have calculated the stochastic background of grav- 
itational waves predicted at second-order due to primor- 
dial density perturbations and scalar metric perturba- 
tions at first order in the early (radiation-dominated) 
universe. In particular we have evaluated the power spec- 
trum of gravitational waves produced first, by a delta- 
function power spectrum, representing an infinitely sharp 
feature in the primordial spectrum at a single wavelength, 
and secondly by a power-law spectrum, compatible with 
direct observations of the primordial density perturba- 
tions on much larger (CMB) scales. It is interesting to 
compare these predictions with the sensitivities of cur- 
rent or planned gravitational wave experiments. 

The resulting power spectrum of gravitational waves 
at matter-radiation equality, r\ eqi can be related to their 
energy density (per logarithmic interval) today by [23j 



£l-y(r]) 47T 
£l<y(r)eq) 3 



(35) 



using the fact that the energy densities of GW and ra- 
diation evolve at the same rate in the matter era. We 
shall consider first the case of the power-law spectrum of 
primordial density perturbations, for which we obtain 



Oy(ri) 124416 a4 



2(n s 



kc 



MB 



T ns (krj eq ) (36) 



We have assumed that all modes are outside the Hubble 
radius at the start of their evolution, io<l 

The latest WMAP data gives /S^(kcMB — 



0.002Mpc _1 ) w 2.36 x 10" 
5 x 10~ 5 today, we have 



GW 



(k,rf) « 2.23 x 1(T 18 



Combining this with fi. 



2(n, — l) 



CMB 



The mode which enters the Hubble radius at matter 
radiation equality defines rj eq — l/k eq where k eq — 
0.009Mpc _1 , corresponding to a GW frequency today of 
f eq = ck eq /2ir sa 1.4 x 10~ 17 Hz. Our pivot value cor- 
responds to frequency Jcmb ~ 3.1 x 10~ 18 Hz. For the 
power-law case we then have 



n GW « 1.86x10 



-20+34(n s -l) 



3.2 



f_ 
'Hz 



where T ns is evaluated at x = oo, and may be esti- 
mated from Fig. [6] Therefore, for a red spectrum with 
n s = 0.95, we find fl G w m 3.2 x 10~ 22 (//Hz) -01 , while 
for a blue spectrum with n s = 1.1 we have flow ~ 
6.1 x 10 -17 (//Hz) ' 2 . In principle, therefore, proposed 



detectors such as LISA 24] or DECIGO [2J, [2j| could 
be used to place limits on the scalar tilt, independently 
of other observations. Of course, we have only consid- 
ered constant tilt here; running of the spectral index will 
change these predictions considerably. Indeed, as the pri- 
mordial power spectrum is directly constrained only at 
very large scales, GWs provide a novel way to probe the 
primordial power spectrum at wavelengths some twenty 
orders of magnitude smaller than CMB measurements 
allow. 

For the case of excess power in a single mode, we may 
define Slew analogously with Eq. (|35p . but near the spike 
Vh is not smooth, so flow is n °t necessarily the actual 
energy density of GW. However, we will use it to give a 
rough indication of power for comparison against quoted 
detector sensitivities. The peak of the spike in the power 
spectrum occurs at f pea k = 2/i„/\/3, which gives 



^Gw(fpeak) « 6.0 X 10~ U A* 1 + 0.09 log 



' 10 lGeV; ' 

We have written this in terms of the temperature at 
which the input wavenumber enters the Hubble radius, 
using [13 f m = ck m /2n » HT 6 (T ent /1 GeV) Hz for 
/ > 10~ 6 Hz, which is the range for any detector of GW. 
We note that the typical resolution of a detector in time 
T is A/ ~ 1/T [23]; so, for a one year observation we have 
A/ ~ 10" 8 Hz. For current detectors such as GEO [H], 
LIGO [3, TAMA [iH or VIRGO 30] whose optimal 
frequency is ~ 100Hz, excess power at a single wavenum- 
ber in the primordial power spectrum would correspond 
to a horizon entry temperature of T ent « 10 8 GeV, im- 
plying that FIgw ~ 10 _16 yl 4 . While this is nominally 
seven orders of magnitude below Advanced LIGO's sen- 
sitivity, Q,gw ~ 10~ 9 , it implies that Advanced LIGO 
could, in principle, pick up modes with A ~ 100, i.e., one 
hundred times the amplitude observed on large scales. 
For detectors such as BBO/DECIGO [Hil], whose op- 
timal frequency is ~ 0.1Hz, their improved sensitivity 
{SIgw ~ 10~ 17 — 10 -15 ) could constrain excess power 
with A ~ 1 at a single wavenumber at a Hubble entry 
temperature ~ 10 5 GeV. 

In summary, we have presented the second-order evo- 
lution equation for tensor modes driven by quadratic 
scalar terms, and have solved this to obtain the resulting 
power spectrum of the gravitational waves. By consid- 
ering density fluctuations with excess power on a given 
scale, as well as a nearly scale-invariant power spectrum, 
we have shown how future gravitational wave detectors 
can place constraints on the primordial density fluctua- 
tions at much smaller scales than can be probed by ob- 
serving the CMB and large-scale structure. 
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APPENDIX A: CO-EFFICIENT MATRICES IN P h 

The coefficients which appear in Eq. (|25p and elsewhere are presented here. First, we define the column matricies 

°=i o i' a = \ +; I- b= i i= c= uj i- (A1) 



Then we may write: 



and, 






^ 1 -^ a + c ) 4,-1(0 



(A2) 



M i ( *f u 2 vl - ±§uv 2 a ^(ua-vl) \ ( 

m nm ~ yx_ u 2 v 2 a o - ± (3u 2 + v 2 )a + \uvl § a ) ' [ ' 

M 2 (k^ 2 b -±{?>u 2 + v 2 )b+\uvc \b\ 

nm \ -^u 2 vc- ^§uv 2 b &(- ub + vc ) o / ' { ' 

M 3 _ ( >§u 2 vl - >§uv 2 a ^(ua-vl) o\ ( 

M 4 _ ( m u2v2b -\{u 2 + v 2 )b+\uvc lb 



